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We propose a simple, general, and accurate formula for analyzing the tunneling between 
classical configurations of a non-planar molecule in a gas medium, as a function of the 
thermodynamic parameters of the gas. We apply it to two interesting cases: i) The shift to 
zero frequency of the inversion line of ammonia, upon increase of the pressure of the gas; ii) 
The destruction of the coherent tunneling of D2S2 molecules in a gas of He. In both cases, 
we compare our analysis with previous theoretical and experimental results. 
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I. INTRODUCTION. 

Quantum tunneling has always been one of the most fascinating — and counterintuitive — 
manifestation of quantum phenomena in nature [U |2]. Among the most well-known examples 
of this kind, are non-planar molecules, in which an atom or a group of atoms oscillates between 
the two wells of the potential energy surface, also at low temperatures [H [3l H]. There has been 
a great amount of interest in the study of coherent tunneling in molecular systems, and how and 
under which physical conditions this coherence may be suppressed due to environmental interac- 
tions [5H1S]. 

The relevant property of the molecule can be effectively modeled by a particle of mass M 
moving in a double- well potential V{q) with two minima at q = ±(70/2; here g is a generalized 
inversion coordinate. The minima, associated with the two localized (say, chiral) states of the 
molecule, are separated by a barrier Vq. If we denote the small-amplitude vibration in either well 
by ujQ, then in the limit Vq ^ hojo ^ ksT, (where T is the temperature of the bath and A;^ is 
Boltzmann constant), the state of the system is effectively confined in the two-dimensional Hilbert 
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space spanned by the two lowest eigenstates. 

As custom, we denote by |+) and |— ) the ground and first excited state of the relevant part 
of the Hamiltonian, separated by the energy difference AE' = hujx- They are delocalized states, 
since they must be also eigenstates of the parity operator. The localized — or chiral — states, are 
the combinations: 

|L> = ^[|+> + |->], |R> = ^[l+>-|->]- (1) 

In this two-dimensional approximation, the Hamiltonian takes the form: H = —uJxd'x/2, where ax 
has I lb) as eigenstates. The position operator q reduces to {qo/2)az, which has the localized states 
|L), |R) as eigenstates. 



II. THE ROLE OF DECOHERENCE. 



A very interesting problem is the behavior of a non-planar molecule in a gas, and how the 
tunneling between different classical configurations depends on the thermodynamic variables of the 
gas. Several approaches have been developed, for analyzing the system-bath interaction. One of 
the best understood ones is decoherence [2l[T71[l8], and in particular collisional decoherence [QtfTOl- 
121] . which is ubiquitous in quantum physics. The application of decoherence to the study of the 
coherent tunneling of a molecule in a gas has been first done by [5l |6]. Different alternative 
approaches to the problem include the mean field theory [TO], or various modeling of the dephasing 
process induced by intermolecular interactions [11], by photons [8J, or by phonons |6]. However, 
decoherence seems to cover the widest ranges of physical situations. 

When analyzing decoherence phenoma, the starting point is the Lindblad-type master equa- 
tion [22J, which in the two-dimensional approximation reads: 

dp ioJx A _ 

where p is the density matrix of the relevant molecule, and A the decoherence rate. The task 
of (collisional) decoherence is to give a microscopic expression for A. In this article we propose 
a general master formula for computing the decoherence rate A for non-planar molecules. The 
advantage of our formula — differently from previous approaches — is that it can be easily applied 
to many different molecules in a large variety of background gases. As we will show, even in its 
simplest form, it gives results which are in excellent agreement with known data. 
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A. A master formula for A. 

The key, and simple, observation is that energy eigenstates are superpositions in space of one 
or more atoms of the molecule, and for this reason they decohere upon interaction with an envi- 
ronment. Therefore we may apply the theory of quantum Brownian motion |20] to such atoms, in 
order to compute the loss of coherence due to the scattering with the bath particles. In the dilute 
gas limit, and when the scattering is recoil- free, the decoherence master equation takes the form 

HMD: 

= -^(x|[^,p]|x') -F(x-x') (xIpIx') (3) 

with H the Hamiltonian of the isolated Brownian particle, |x) , |x') its position eigenstates, and: 

P 

F (x) = n dpv {p) — a {p, x) (4) 
Jo 1^ 

where n is the density of bath particles, m the reduced mass of bath particles, p the bath particle 
momentum, i/ {p) the distribution of such momenta, and: 

a(p,x) = |^^(l-e-(--"^)-/^^)|/(p,^)|^ (5) 

is the decoherence cross section; / (p, 9) is the scattering amplitude, and 9 the scattering angle 
between directions fii and n2. 

In the limit where the two-dimensional approximation is valid, the spatial coherence is between 
the two minima of the double-well potential, which are separated by a distance |x — x'| = qq. We 
assume that, for all important bath particles' momenta the decoherence cross section a (p, qq) 
does not change appreciably. We can then write: 

Tip 

A = F{qo) ~ —cr{p,qo), (6) 
m 

where p = \/2mkBT is the most probable momentum, for an ideal gas at equilibrium at temperature 
T. For isotropic scattering, we have [20] : 

a {p, go) = 2^ d (cos 9)[l-^-^^\f{p,9)\^ (7) 

with e = {2pqQ/h)s\u{9/2). Eq. can be easily computed in two limiting cases: the low 
temperature limit pq^ <^ h, i.e. T <^ h? /2mkBqQ', and the high temperature limit pqo S> h, 
i.e. T > h^/2mkBql. 

For values of ojq in the infra-red range (10^3-10^^ Hz), which is the typical range of molecular 
vibrations, go of the order of few Angstrom, and m of 2-30gr/mol for masses of typical background 
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gases (e.g., from H2 to Air), one has: Huo/kB ~ 500-5000 K and It? /2mkBqQ ~ 1-30 K. This 
means that both the two dimensional approximation and the high temperature limit are valid 
in most interesting physical situations. In this case, the variable is large, sin 0/0 ~ 0, thus: 
c {Pi Qo) = cttot (P) = 2Tr J d (cos0) \ f {p, 9)\'^, which is the total cross-section. 

If we consider the scattering interaction as given by a central potential and use the 
phase-shift method 121! , the scattering amplitude / {p, 6) can be expressed as: / (p, 9) = 
{h/p) X^^o (^^ + ^) sinSi Pi {cos0), with (cos 6*) the Legendre polynomials, and 6£ the phase 
shifts. The phase shifts have suitable integral representations depending on the scattering potential. 
Collecting the above results, we can write: 

where we used n = P/ksT. 

Eq. ([s]) and, with it, decoherence theory, is very powerful: it tells that the effect of a gas medium 
on the tunneling properties of a molecule can be reduced to the two-body interaction (scattering) 
between the relevant molecule and a single gas molecule. The gas as a whole enters only though 
its thermodynamics parameters. Eq. ([s]) is also very general, because can be immediately applied 
to any situation, where a non-planar molecule is immersed in a gas medium (within the limits 
of validity of the general approach, previoulsy discussed). Eq. ([s]) is also independent from the 
internal- rotational states of the non-planar molecule. 

Computing the phase shifts 61 can prove quite hard. One can provide a handy expression in 
the hard-sphere limit: sii? 6e = jj {pa/h) /[jf (pa/h) + nj (pa/h)] with a the distance of closest 
approach of the centers of the two colliding spheres, ji (x) the spherical Bessel functions and n£ (x) 
the spherical Neumann functions j23|. [M] . The hard-sphere limit is the zero-th order approximation 
one can make. However it turns out to be very good, because intermolecular interactions typically 
have a finite effective range, and the main contribution to A comes essentially from scattering in 
this effective region. Considering the effective range as that of a hard-sphere, we can find a very 
simple and general expression for A. In fact, as we will show in the following examples, the hard- 
sphere limit provides a very good description of decoherence effects in the tunneling properties 
of a molecule. It also offers a geometric, and very intuitive, way to understand the decoherence 
mechanism, reducing it to the bouncing of two spheres. In this way, one can understand why, in the 
calculations, one has to include a steric factor 7, which takes into account only those directions of 
approach of a gas molecule, which can contribute to decoherence: an atom attached to a molecule 
can be "hit" by a bath particle only from the directions which are not hindered by the rest of the 
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molecule. We now apply this setup to two interesting physical situations. 



III. SHIFT IN INVERSION-LINE FREQUENCY OF AMMONIA GAS. 

In [25], the inversion- line frequency of ammonia has been experimentally measured as a function 
of the gas pressure P. The frequency is found to decrease from 0.78 cm^^ at P = 0, to zero above 
P = 2 atm (See Fig. [T]). This phenomenon was first explained qualitatively in \26\ I27j ; to our 
knowledge, the first complete theoretical analysis is provided in [lOj, based on the mean field 
theory. The inversion frequency turns out to be: 



i^x = yi - (9) 

where P^r = 1.695 atm is the critical pressure derived theoretically. The theoretical formula, re- 
ported in Fig. [T| is in very good agreement with the experimental data. 

The mean field theory provides a 'static' approach to the problem. Decoherence theory provides 
a 'dynamical' approach, and its predictions partly differ, as we now show. This comes through 
the explicit solution of Eq. By writing p = (1/2) (/ + Pxd'x + Py^y + PzS'z), all information 
is contained in the Bloch vector P = {Px,Py,Pz)- The experimental results obtained by cavity- 
resonator technique are connected to the relative occupation probability Pz (t) [25J, whose time 
evolution, in the under-damped case, reads [28] : 

Pz{t) = —^[2PziO)uJx COS U!xt (10) 
ZUJx 

+ {Pz{0)X + 2Py {0)cOx) sin Cdxt], 
where the characteristic inversion frequency, Ox, is: 



px2 



Ox = Ux^l-^-J , (11) 

and Per = 2Pu;x/X- Since the experiment in has been performed at room temperature, the high 
temperature limit holds and, according to Eq. ([s]), we can write the critical pressure as follows: 

P. = V^M^, (12) 



As expected, Eq. (11) predicts a decrease of the inversion frequency, due to the interaction of each 



molecule with the rest of the gas. Moreover, Eq. (12) provides a microscopic expression for P„, 
without free parameters. 
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FIG. 1: Inversion frequency of ammonia, as a function of the gas pressure P. The dots are experimental 
data taken from [3S]; the solid line refers to Eq. (11), with Per = 1.05 atm; the dashed line refers to Eq. ^ 
of [TU], with Per = 1.695 atm. In both cases, ujx (frequency at P = 0) has been taken equal to 0.78 cm~-^, as 
reported in |25j . 



A good estimate for the total cross section (Ttot(p) can be provided in the hard-sphere limit, 
previously discussed. The hard-sphere radius is equal to twice the effective hard-sphere radius 
of ammonia (a = 4.38 A) [29]. In computing the total cross section, the relevant contribution 
comes from those terms with i < pa/h ~ 22 [23l [23]; this gives a^oT = 481a|^^^ (where Oeohr 
is Bohr radius). In the calculations, we have chosen the masses as follows. According to the 
double-well description, the mass M of the Brownian particle corresponds to the reduced mass 
M = g^"^^^ pOj with mil niass of Hydrogen and mass of Nitrogen. The mass m entering 
the calculation of the total cross section and of the decoherence rate refers to the reduced mass of 
the Brownian and bath particle (m* in Eq. (2.3) of [20]): m = ^^^^^^ , with rriNHg the mass of 
ammonia. We have to multiply the cross section by a steric factor 7 ~ 0.85, since not all collisions 
of a bath particle with the Brownian particle result to decoherence (see Fig. 2(a)). By taking 



LJ^ = 0.78 cm-^ and T = 300K [25], Eq. (|12|) gives: Per = 1.05 atm. 

As shown in Fig. [T| Eq. (11) matches well the experimental data at low pressures, but clearly 
fails at high pressures. This does not come as a surprise, since the microscopic model used to 
compute A is valid only in the diluted gas case. As explained in [3], in the case of ammonia, the 
region P > 1 atm, where many-body collisions become important, is far from the dilute gas limit. 

It is important to note that Eq. (11) and Eq. ^ predict a different behavior at very low pres- 



sures. While, according to Eq. (11), the inversion frequency Cox approaches cOx with an horizontal 
slope for P — 7- 0, according to Eq. ([9]) the slope is negative. The experimental data seem to suggest 
that the slope gradually becomes horizontal, confirming our prediction. 
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(a) 



(b) 



FIG. 2: (a) Collision between an NH3 Brownian molecule (left) and a NH3 bath particle (right); we simplify 
the description by considering the bath particles as spheres. When a bath particle collides with the Brownian 
particle within the solid angle depicted in the picture, there is no significant decoherence effect, because 
the H atoms (who effectively participate to the superposition) are not hit. Accordingly, the geometric cross 
section must be multiplied by a steric factor 7 — {An — il)/47r. This takes into account only those collisions 
which give rise to decoherence. With simple geometric considerations, one can estimate the angle 6 which 
identifies fl as: 9 c^i 1.088. Then Q, — 2t: sinxdx ~ 1.875, and 7 ~ 0.851 (b) Collision of a Brownian 
D2S2 molecule (right) with a bath He particle (left). In this case, only the D atoms participate to the 
superposition. Each D atom is hit by a bath particle when this arrives within the solid angle fl shown in the 
picture. Simple geometric considerations lead to the result: ~ 4.943. The steric factor for each D atom 
then is 7 = 0/47r ~ 0.393. 



There is also a strong theoretical reason why it should be so. Eq. ( 11 ) is a direct consequence of 



Eq. ([2]), i.e. of a Lindblad-type master equation, which is the result of very general mathematical 
requirements, such as the quantum dynamical semigroup (QDS) structure, and complete positiv- 
ity [22]. Thus, it would be interesting to perform the experiment again, to assess the low pressure 
behavior of the inversion-line frequency of ammonia. 



IV. EFFECT OF DECOHERENCE ON TUNNELING DYNAMICS OF D2S2 

MOLECULE. 

In [13] the decoherence rate A of a bath of Helium (He) on a D2S2 molecule is studied with heavy 
numerical machinery. Actually, its value is not explicitly given, but can be estimated from Eq. (2) 
of the paper. We assume that the decoherence cross section ijtot (notation of [13j) is constant over 
significant bath particles' momenta and, according to Fig. (1) of the paper, equal to lOSal^,^^ for 
the considered temperature T = 300 K. Then Eq. (2) of the paper gives A = 215.5 Hz. 

We now apply our framework to give an independent estimate for A, using Eq. ([s]). The Brownian 
particle has the reduced mass M = "^^."^^ with mo the mass of Sulfur and m-n the mass of 
Deuterium. For m we use m = with rrin^ the mass of He. In the hard-sphere limit, it is 
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reasonable to consider a as the sum of the hard-sphere radius of He (1.085 A) and the van der 
Waals radius of D (1.2 A) [29]. Under the conditions of P = 1.6 x 10"^ Pa and T = 300 K used 
for the experiment proposed in [13| . the high temperature hmit holds. In computing the total 
cross section, the relevant contribution comes from those terms with i < pa/h ~ 9 [23l I24| . One 
obtains crTOT(p) = 139. 4a|^^^. This result has to be multiplied by a factor 2, since there are two 
D atoms superimposed for each molecule, and also by a steric factor 7 ~ 0.393, which takes into 
account that each D atom does not collide with bath particles coming from all directions, due to 



the presence of the sulfur atoms (see Fig 2(b)). This gives a-ro-rip) = 109.7a|^jj^ and, according to 
Eq. ([s]), A = 229.4 Hz. Our result is in very good agreement with the one predicted by [Iffj . 



V. CONCLUSIONS. 



Our approach to the tunneling properties of non-planar molecules in a gas offers a simple and 
general formula for computing the decoherence rate, and its dependence on the thermodynamic 
variables of the gas. It is flexible and applies to all situations where a molecule is immersed in a 
gas, which is the most common case, when tunneling properties are under study. The only three 
limitations are: the diluted-gas limit, the two-dimensional approximation and the high-temperature 
limit, which are satisfied in most physical conditions. The hard sphere limit represents the grossest 
type of approximation, but as we have seen its predictions are in excellent agreement with previous 
analyses. And it has the advantage of reducing a complex many-body problem to a simple geometric 
calculation. This opens the way to easy estimates of decoherence effects also for very complicated 
molecules in highly structure gases. 
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